Astronomy & Astrophysics manuscript no. mixpns 


©ESO 2013 


February 26, 2013 





Structure of the hadron-quark mixed phase in protoneutron stars 

H. Chen, 1 - 2 G. F. Burgio, 1 H.-J. Schulze, 1 - 3 and N. Yasutake 4 



1 INFN Sezione di Catania, Dipartimento di Fisica, Universita di Catania, Via Santa Sofia 64, 95123 Catania, Italy 

2 Physics Department, China University of Geoscience, Wuhan 430074, China 

3 



Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 
4 Department of Physics, Chiba Institute of Technology, 2-1-1 Shibazono, Narashino, Chiba 275-0023, 

February 26, 2013 



Japan 



o 

(N 

in 

(N 



ABSTRACT 

We study the hadron-quark phase transition in the interior of hot protoneutron stars, combining the Brueckner-Hartree-Fock approach 
for hadronic matter with the MIT bag model or the Dyson-Schwinger model for quark matter. We examine the structure of the mixed 
phase constructed according to different prescriptions for the phase transition, and the resulting consequences for stellar properties. 
We find important effects for the internal composition, but only very small influence on the global stellar properties. 
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T71. Introduction 
PS 

^\ A protoneutron star (PNS) is formed after the gravitational col- 
JH lapse of the core of a massive star (M > 8M Q ), exploding in a 
P-itype-II supernova (Shapiro & Teukolsky 1983; Bethe 1990). Al- 
q though the explosion mechanism is still not fully explained (Bur- 
$_i rows 2012), some general features can be considered as robust. 
^ In fact, just after the core bounce, the PNS is very hot and lepton- 
C3 rich, and neutrinos are trapped for a few seconds. The following 
l— ~~ 1 evolution of the PNS is dominated by neutrino diffusion, caus- 
! ing deleptonization and subsequently cooling. Ultimately, the 
neutron star (NS) achieves thermal equilibrium and stabilizes at 
practically zero temperature without trapped neutrinos. 
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The theoretical description of the formation of a PNS re- 
quires an accurate treatment of the microphysics of the collaps- 
• ing matter, in particular of neutrino transport and related pro- 
cesses (Burrows & Lattimer 1986; Prakash et al. 1997; Burrows 



o 



2012). Moreover the violent dynamical processes occurring in 
the contracting-exploding star need to be treated in the frame- 
. . work of general relativity [see Ott (2009) for a recent review]. 
J> The physical processes which contribute to the subsequent PNS 
K^J evolution, such as nuclear and weak interactions and energy and 
l^j lepton number transport by neutrino diffusion, are very difficult 
to include in dynamical simulations. Thus, most simulations of 
the gravitational core collapse to a PNS end shortly after the core 
bounce and the launch of the supernova explosion, typically af- 
ter a few hundreds of milliseconds, and only a few dynamical 
simulations extend to the first minute of the PNS life (Pons et al. 
1999; Fischer & Mueller 2009). 

During the evolution of a PNS into a NS, a hadron-quark 
(HQ) phase transition could take place in the central region of 
the star (Prakash et al. 1995; Lugones & Benvenuto 1998; Pons 
et al. 1999, 2001; Steiner et al. 2000; Epsztein Grynberg et al. 
2000; Nicotra et al. 2006a; Yasutake et al. 201 1), and this would 
alter substantially the composition of the core. In fact the heav- 
iest NS, close to the maximum mass (about two solar masses), 
are characterized by central baryon densities larger than 1 /fm 3 , 



as predicted by calculations based on a microscopic nucleonic 
equation of state (EOS). 

The study of hybrid stars is also important from another 
point of view: Purely nucleonic EOS are able to accommodate 
fairly large (P)NS maximum masses (Baldo et al. 1997; Akmal 
et al. 1998; Glendenning 2000; Zhou et al. 2004; Li & Schulze 
2008), but the appearance of hyperons in beta-stable matter 
could strongly reduce this value (Glendenning 2000; Schulze 
et al. 2006; Li & Schulze 2008; Carroll et al. 2009; Dapo et al. 
2010; Schulze & Rijken 2011). In this case the presence of 
non-baryonic, i.e., "quark" matter would be a possible manner 
to stiffen the EOS and reach larger NS masses (Burgio et al. 
2002; Maieron et al. 2004; Kurkela et al. 2010; Weissenborn 
et al. 201 1). Heavy NS thus would be hybrid quark stars. 

In previous articles (Nicotra et al. 2006b; Burgio & Schulze 
2009, 2010; Burgio et al. 2011a) we have studied static prop- 
erties of PNS using a finite-temperature hadronic EOS includ- 
ing also hyperons (Burgio et al. 2011b) derived within the 
Brueckner-Bethe-Goldstone theory of nuclear matter (Baldo 
1999). An eventual HQ phase transition was modeled within 
an extended MIT bag model (Nicotra et al. 2006a; Yasutake 
et al. 2011) or a more sophisticated quark model, the Dyson- 
Schwinger model (DSM) (Roberts & Williams 1994; Roberts & 
Schmidt 2000; Alkofer & von Smekal 2001; Roberts et al. 2007; 
Chen et al. 2011,2012). 

The purpose of the present work is to complement our pre- 
vious articles by studying details of the HQ phase transition oc- 
curing in hybrid stars and their implications for the structure of 
a PNS, in particular the question whether global (P)NS observ- 
ables are sensitive to and thus may reveal information on the 
internal stellar structure. 

In Sec. II we briefly sketch the theoretical approaches which 
we use for modeling the hadron and the quark phases, and in 
Sec. Ill we describe the corresponding pure phases. The struc- 
ture of the mixed phase is discussed in Sec. IV, and the results 
are illustrated in Sec. V. Finally, in Sec. VI we summarize our 
conclusions. 
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2. Equations of state 



servation: 



The EOSs for hadronic matter (HM) and quark matter (QM) that 
we use in this work, have been amply discussed in previous pub- 
lications (Burgio et al. 2002; Maieron et al. 2004; Nicotra et al. 
2006a,b; Burgio & Schulze 2009, 2010; Burgio et al. 2011b; 
Chen et al. 2011, 2012), where all necessary details can be 
found. Our hadronic EOS is obtained from Brueckner-Hartree- 
Fock (BHF) calculations of (hyper)nuclear matter (Schulze et al. 
1998; Baldo et al. 1998, 2000) based on realistic potentials [the 
Argonne Vis nucleon-nucleon (Wiringa et al. 1995) and the Ni- 
jmegen NSC89 nucleon-hyperon (Maessen et al. 1989) in this 
case] supplemented by nucleonic Urbana U1X three-body forces 
(Carlson et al. 1983; Schiavilla et al. 1986; Pudlineret al. 1997), 
and extended to finite temperature (Burgio et al. 201 lb). We em- 
ploy two different representative models for QM, an extended 
MIT bag model [the model with a density-dependent bag con- 
stant of Burgio et al. (2002); Maieron et al. (2004); Nicotra et al. 
(2006a)] and a Dyson-Schwinger model [the model DS4 of Chen 
et al. (2011, 2012)], which yield in fact quite different internal 
structures of hybrid stars. 

Those theoretical calculations provide the free energy den- 
sity of the bulk system (pure HM or QM) as a function of the rel- 
evant partial number densities n, and the temperature, /({«,}, T), 
from which all thermodynamic quantities of interest can be com- 
puted, namely, the chemical potentials yu„ pressure p, entropy 
density s, and internal energy density s read as 



n B 



Hi 
P 

s 
s 



= n 



drii ' 

2 d(f/n B ) 



dT ' 
f + Ts, 
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(3) 
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where n B is the total baryon number density. These quanti- 
ties allow to determine the stellar matter composition and the 
EOS, which is the fundamental input for solving the Tolman- 
Oppenheimer-Volkoff equations of (P)NS structure. 



3. Pure phases 

In neutrino-trapped beta-stable (hyper)nuclear or 
quark matter the chemical potential //, of any particle 
i = n,p, A, IT, u, d, s, e,fi, v e , v^, . . . is uniquely determined 
by the conserved quantities baryon number B„ electric charge 
C„ and weak charges (lepton numbers) L®, I - e,p. with 
the corresponding set of independent chemical potentials 



[ii = Bijj. B + Que + Lf } /j. L <e) + L^Vlm 



r(M), 



(5) 



In this work we neglect muons and muon neutrinos due to their 
low fractions and negligible impact on global stellar properties, 
hence use simply v = v e , L = L (e \ The relations between chem- 
ical potentials and partial densities for hadrons and quarks are 
given by the microscopic models mentioned before, while lep- 
tons are treated as free fermions. With such relations, the bulk 
system in each phase can be solved for a given baryon density, 
imposing the charge neutrality condition and lepton number con- 



o = ^md, 

i 

Y e n B = Yj 1 ^ 



(6) 
(7) 
(8) 



When the neutrinos v e are untrapped, the lepton number is not 
conserved any more, the density and the chemical potential of v e 
vanish, and the above equations simplify accordingly. 



4. Mixed phase constructions 

We are interested in the HQ phase transition in PNS and consider 
therefore the usual oversimplified standard conditions, namely 
trapped hot matter with a fixed lepton fraction Y e = (n e + 
n v )/n B = 0.4 and either an isentropic, S /A = 2, or an isother- 
mal, T = 40 MeV, temperature profile. One could consider more 
realistic profiles (Burgio et al. 201 la), but we focus in this work 
on the difference between phase transition constructions. 

A fully microscopic treatment of the HQ mixed phase involv- 
ing finite-size (pasta) structures can only be performed numeri- 
cally (Tatsumi et al. 2003; Endo et al. 2005, 2006; Maruyama 
et al. 2007; Yasutake et al. 2009, 2012a,b). One introduces 
Coulomb energies and surface energies via a HQ surface tension 
and then minimizes the (free) energy of a Wigner-Seitz (WS) 
cell, allowing for different geometrical structures of the quark 
phase embedded in the hadron phase and vice versa. The output 
are the optimal size and geometry of the cell, as well as the lo- 
cal distributions of the individual particle species, and also the 
Coulomb field inside the cell. Some illustrative examples can be 
found in the given references. 

This is a very time-consuming and not very transparent nu- 
merical procedure. It is therefore convenient to search for re- 
liable approximations to this procedure, and in this article we 
compare two prescriptions corresponding to two limiting cases 
of the full numerical procedure, that are termed global charge 
neutral (GCN) and local charge neutral (LCN) mixed phase. 

The first procedure is well known as Bulk Gibbs or Glen- 
denning construction (Glendenning 1992, 2001) from the zero- 
temperature case and corresponds to a "small" WS cell [com- 
pared to the electromagnetic Debye screening length, which is 
about 5-10 fm (Heiselberg et al. 1993; Heiselberg 1993; Glen- 
denning & Pei 1995; Christiansen & Glendenning 1997; Takat- 
suka et al. 2006)] caused by a "small" HQ surface tension. In 
this case the electromagnetic potential is practically constant 
throughout the cell, and an electric field does not exist. Con- 
sequently the electron density is also constant, while the hadron 
and quark densities and their electric charges are different in or- 
der to fulfill the conditions of pressure and baryon chemical po- 
tential equality at the HQ interface. In the case of neutrino trap- 
ping, the neutrino densities have also to be equal in both phases, 



(9) 



which together with the equal electron densities implies equal 
lepton densities tii = n e +n v (but not lepton fractions Y e = niln B ) 
in both phases. Altogether we have therefore the equality of the 
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Fig. 1. Relative populations x, = iij/n B for different stellar composi- 
tions: trapped matter with S /A = 2 (upper panels) or T = 40 MeV (cen- 
tral panels), and untrapped matter with T = (lower panels). The GCN 
(left panels) or LCN (right panels) construction for the mixed phase is 
employed together with the MIT bag model for the quark phase. 



intensive thermodynamical quantities in both phases: 
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which together with the general rule Eq. (5) determines the com- 
position of the system for given overall baryon density n B , van- 
ishing electric charge, fixed lepton fraction Y e in the trapped 
case, and eventually a prescribed entropy profile S IA{n£): 



(\-X)n H B +Xn Q B 

(1 - X )n H L + X n Q L 
(\-X)s H + X s Q 



= n B , 

= 0, 

= n B Y e , 

= n B S/A. 



(15) 
(16) 
(17) 
(18) 



where x is the volume fraction occupied by the quark phase and 
the last equation determines the local temperature. 

The opposite limiting case (LCN) corresponds to a WS cell 
that is large relative to the electromagnetic Debye screening 
length, and a large surface tension. In this situation the elec- 
tric charges are well screened inside the cell and both QM and 
HM are locally charge neutral nearly everywhere, 

n£ = «2 = 0, (19) 

except on a small boundary layer near the HQ interface, where 
a positively charged layer of HM and a negatively charged one 
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Fig. 2. Same as Fig. 1, but with the Dyson-Schwinger model for the 
quark phase. 



of QM are present and create a strong but very localized elec- 
tric field (Voskresensky 2002; Voskresensky et al. 2003). Con- 
sequently there occurs a sharp rise 6^c of the Coulomb potential 
at the HQ interface and Eq. (1 1) is modified to 



(20) 



such that for example the electron density is now different in 
hadron and quark phases. 

In beta-stable untrapped matter this situation corresponds ex- 
actly to the usual Maxwell construction, joining two charge- 
neutral phases by equality of pressure and baryon chemical po- 
tential, that is often employed for simplicity. Including neutrino 
trapping with microscopic finite-size structures requires always 
homogeneous neutrino densities, Eq. (9), and therefore, due to 
the unequal electron densities, in this case the trapping condi- 
tion becomes a global one, as expressed by Eq. (17). Due to the 
additional degree of freedom represented by the neutrino density, 
then the LCN construction is realized in the PNS as an extended 
mixed phase involving a HQ coexistence region with a continu- 
ously varying pressure (Hempel et al. 2009; Pagliara et al. 2009, 
2010; Yasutake et al. 2012b). 

We have explained that the GCN and LCN constructions are 
in fact idealized scenarios that correspond to two opposite ex- 
tremes of the microscopic treatment. It has been pointed out in 
Yasutake et al. (2012a,b) that actually the LCN construction is 
closer to the full microscopic treatment of finite-size effects than 
the GCN, and it is therefore of interest to compare the predic- 
tions of the two constructions for the internal composition and 
other properties of PNS, which we will do now. 
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5. Results 

5.1. Internal composition 

The relative particle populations are shown as a function of the 
baryon density in Fig. 1 for the bag model and in Fig. 2 for the 
DSM, for trapped matter with Y e = 0.4 and i) entropy per baryon 
S I A — 2 (upper panels), ii) temperature T — 40 Me V (middle 
panels), and iii) untrapped and cold neutron star matter (lower 
panels). The GCN (left panels) and LCN (right panels) calcula- 
tions are compared. 

There are big differences between the MIT and DSM regard- 
ing the HQ mixed phase that have been pointed out in Chen et al. 
(201 1, 2012): With the MIT model the HQ phase transition starts 
at fairly low baryon density and a pure quark phase is reached at 
not too large density, whereas with the DSM the onset of the 
mixed phase occurs at higher density and the system remains in 
the mixed phase even at very large density. Furthermore, hyper- 
ons are allowed with the MIT model, where they might appear 
only in small fractions at low density and are replaced by strange 
QM at higher density (this can be seen in the central panels of 
Fig. 1), whereas they prevent any transition to QM with the DSM 
and have to be excluded by hand in that case. Therefore, these 
two very different quark models might be good candidates to 
reveal important differences between the phase transition con- 
structions that we are examining. 

In fact, comparing now the results obtained with both pre- 
scriptions (left and right panels), we observe behavior in line 
with the general properties mentioned before. We remind that 
for cold NS matter (bottom panels) the LCN corresponds to the 
usual Maxwell construction and the GCN to the bulk Gibbs con- 
struction; and it is well known that the density range of the mixed 
phase with the Gibbs construction is wider than the one with the 
Maxwell construction, which can be seen in the plots. This be- 
havior remains also in trapped hot matter, where the GCN spans 
always a wider density range than the LCN (with the MIT) or 
begins at lower density (with the DSM). In all cases the trapping 
condition shifts the onset of the mixed phase to slightly higher 
density. 

The differences between the LCN and GCN constructions 
are fairly small for the MIT model, but significant for the DSM: 
Here the GCN (Bulk Gibbs) mixed phase occurs in a much wider 
density interval than the LCN (Maxwell) one. Apart from the 
Maxwell construction for cold NSs, the matter remains in the 
mixed phase and pure QM is never reached. This variance is due 
to the qualitatively different density dependence of the effective 
bag constant in the MIT and DSM, see Chen et al. (201 1, 2012). 

5.2. Equation of state 

In Fig. 3 the EOS p(ns) is displayed for the different stellar con- 
figurations, quark models, and mixed phase constructions as be- 
fore. For comparison also the pure phases (nucleons only, nucle- 
ons+hyperons, quarks) are shown. 

We observe only minor differences between GCN and LCN 
in the hot and trapped matter, even for the DSM, where the par- 
ticle fractions are quite different in both cases. Stronger differ- 
ences between GCN and LCN appear in the cold case, where a 
plateau in the pressure shows up for the LCN (Maxwell) calcu- 
lation. 

Thus if during the temporal evolution the system would re- 
main in the LCN phase, the extended mixed phase region would 
gradually disappear. However, one should remember that both 
LCN and GCN are highly idealized constructions and the true 
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Fig. 3. The pressure vs. baryon density for the same stellar configura- 
tions as in Fig. 1, and obtained with the MIT (left panels) and the DSM 
(right panels) quark models. 



behavior depends on currently uncertain microphysics like the 
HQ surface tension. 

5.3. Stellar structure 

Once the EOS is known, the stable configurations of a (P)NS can 
be obtained from the well-known hydrostatic equilibrium equa- 
tions of Tolman, Oppenheimer, and Volkov (Shapiro & Teukol- 
sky 1983). In the low-density range, where nucleonic clustering 
sets in, we cannot use the BHF approach, and therefore we join 
(Burgio & Schulze 2010) the BHF EOS to the finite-temperature 
or -entropy EOS of Shen et al. (1998a,b), which is more appro- 
priate at densities below ris < 0.07 frrT 3 , since it does include 
the treatment of finite nuclei. 

Our results for the gravitational mass as a function of the 
central baryon density, for the different stellar configurations and 
using the different EOS introduced previously, are displayed in 
Fig. 4. We observe in NS matter the strong softening effect of 
the hyperons (dash-dotted purple curves) on the purely nucleonic 
configurations (dotted blue curves), which is however strongly 
reduced in trapped matter, because the hyperon concentrations 
remain smaller (Prakash et al. 1997; Burgio et al. 201 lb). 

In the case of the MIT model, the mass-density relations 
obtained with the GCN or LCN constructions (dash-dot-dotted 
green vs short-dashed orange curves) are nearly indistinguish- 
able, apart from the unphysical low-mass region, where the 
Maxwell construction can be recognized in the NS configura- 
tion with LCN. The maximum mass of the hybrid stars decreases 
slightly with respect to both the nucleonic and the hyperonic 
stars in the trapped cases, whereas it increases (decreases) with 
respect to cold hyperon (nucleon) NS. In NS this is due to the 
fact that the hyperon population is suppressed by the onset of 
quarks, whereas in PNS the trapping reduces the hyperon pop- 
ulation. In all cases, with the HQ phase transition, the value of 
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Fig. 4. The gravitational mass vs. central baryon density for different 
stellar configurations and EOS. 



the maximum mass is about 1.5 M Q and thus rather low, as is a 
general feature of the MIT model (Alford & Reddy 2003). 

If the DSM is used, the differences between GCN and LCN 
(solid black vs dashed red curves) are slightly larger, in partic- 
ular for NS, where the LCN leads to unstable configurations at 
the onset of the quark phase. Nevertheless the differences be- 
tween GCN and LCN maximum masses are insignificant in all 
configurations. In this case the phase transition takes place only 
if hyperons are excluded from the hadronic phase, and the value 
of the maximum mass for the hybrid configuration always de- 
creases with respect to the purely nucleonic star. In particular, 
for PNS the maximum mass is about 1 .75 M , while for NS it is 
slightly smaller. 

These values for the maximum mass of hybrid (P)NS depend 
on the choice of the nucleonic three-body force and on the de- 
tails of the DSM. Larger values could eventually be reached with 
different parameter choices (Chen et al. 2011, 2012), in agree- 
ment with the current observational data (Demorest et al. 2010), 
which is in contrast to the case with the MIT model. However, 
the exact value of the maximum mass of a (P)NS is still an open 
problem [see also Romani et al. (2012)] and not the purpose of 
the present article. 



6. Conclusions 

In this article we discussed the occurrence of a hadron-quark 
mixed phase in the interior of hybrid (proto)neutron stars. We 
explained the physical origin and justification of the idealized 
LCN and GCN phase transition constructions, which represent 
two opposite limiting cases of the microscopic treatment of elec- 
tromagnetic finite-size effects, and examined their consequences 
with two very different quark models. 

While indeed the internal composition of hybrid (proto)stars 
turns out to be very different with both constructions, the impact 
on the equation of state and masses is very much reduced, so that 
these global observables could hardly serve as an indication for 
the type of phase transition and thus the internal stellar structure. 

Therefore, for a true understanding of the nature of the mixed 
phase, detailed microscopic investigations of the finite size ef- 
fects and their importance for the stellar microphysics (cooling, 
transport, oscillations, ...) are required. 
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